QCD phase diagram from finite energy sum rules 
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We study the QCD phase diagram at finite temperature and baryon chemical potential by relating 
the behavior of the light-quark condensate to the threshold energy for the onset of perturbative QCD. 
These parameters are connected to the chiral symmetry restoration and the deconfinement phase 
transition, respectively. This relation is obtained in the framework of finite energy QCD sum rules 
at finite temperature and density, with input from Schwinger-Dyson methods to determine the light- 
quark condensate. Results indicate that both critical temperatures are basically the same within 
some 3% accuracy. We also obtain bounds for the position of the critical end point, /i_g c > 300 MeV 
and T c < 185 MeV. 
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I. INTRODUCTION 

In Quantum Chromodynamics (QCD) the strong in- 
teraction among quarks depends on their color charge. 
When quarks are placed in a medium this color charge 
is screened with increasing density. The density can in- 
crease either by raising the temperature, so that collisions 
between quarks produce more quarks and gluons, or by 
compressing the system, thereby increasing the baryon 
density. If the density increases beyond a certain critical 
value one expects that the interactions between quarks 
no longer confine them inside a hadron, so that they are 
free to travel longer distances and deconfmc. This tran- 
sition from a confined to a deconfined phase is usually 
referred to as the deconfinement phase transition. 

A separate phase transition takes place when the 
realization of chiral symmetry shifts from a Nambu- 
Goldstone phase to a Wigner-Weyl phase. In the mass- 
less quark limit this is achieved by the vanishing of the 
quark condensate, or alternatively the pion decay con- 
stant. Qualitatively, one expects these two phase transi- 
tions to take place at approximately the same tempera- 
ture. An outstanding issue is whether this conclusion also 
holds quantitatively. To address this, it has been custom- 
ary to study the behavior of the order parameters of these 
transitions as functions of the temperature T and the 
baryon chemical potential fis, namely the Polyakov loop 
L [2j and quark anti-quark condensate (^'0) in the chiral 
limit, respectively. In the confined phase the former pa- 
rameter either vanishes in the limit of massless quarks, or 
else it is exponentially suppressed for finite quark masses, 
while it is finite in the deconfined phase. The quark con- 
densate is finite in the confined phase, while it vanishes in 
the deconfined phase at high enough temperature, and in 
the limit of massless quarks. For finite quark masses chi- 
ral symmetry is explicitly broken at the Lagrangian level 



and therefore the phase transition is suppressed. This is 
similar to what happens to a ferromagnet in the presence 
of an external magnetic field. In this situation one might 
need to specify to what extent one is still dealing with a 
phase transition. 

At finite T, and ijlb = 0, lattice QCD calculations pro- 
vide a consistent quantitative picture of the above be- 
havior, resulting in similar critical temperatures T c for 
both transitions in the range 170 MeV < T c < 200 MeV, 
for finite quark masses 3 5]. The situation is much less 
clear cut when both T and /jb arc simultaneously non- 
zero. Lattice QCD simulations cannot be used for [Lb 7^ 
because the fermion determinant becomes complex and 
thus standard Monte Carlo methods fail, as the integrand 
is no longer real and positive definite. However, these 
techniques can still be adapted to extract some, though 
not exact, information on the QCD phase diagram for 
Ms 7^ [3|. Therefore, one needs to resort either to 
mathematical constructions to overcome the above lim- 
itation Q, or to model calculations Of particular 
recent interest is the search for a possible critical end 
point that signals the strengthening of the order of the 
transition with increasing /ig, indicating the beginning of 
a true chiral symmetry restoring/deconfining phase tran- 
sition. The results from Monte Carlo simulations and 
model calculations, with and without Polyakov loop, or 
its variants, seem to be in conflict. In fact, the former give 
smaller (larger) values for the end point baryon chemi- 
cal potential (temperature) than the latter. Things be- 
come worse if one uses the imaginary chemical potential 
method, a well established technique for not too large 
values of [Ib- Indeed, this leads to a shrinking and weak- 
ening region of chiral phase transitions with increasing 
Hb, thus suggesting that there is no critical end point 
f° r Ms ^ 500 MeV [Toj j . It has also been pointed out 
that even if the transition weakens with increasing fis. 
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the existence of the critical end point would not be ruled 
out, although it would require a non-monotonic behav- 
ior [ll|. In view of this situation, alternative ways of 
examining the QCD phase diagram are required. 

One possibility is to look at variables that describe 
deconfmement other than the Polyakov loop. A phe- 
nomenological QCD parameter associated with decon- 
fmement was first proposed long ago in [l2j . and it is 
the square energy (so) beyond which the hadronic res- 
onance spectral function becomes smooth and well de- 
scribed by perturbative QCD (PQCD). At T = this 
continuum threshold lies in the range s — 1 — 3 GeV 2 , 
depending on the channel. At finite temperature one 
expects So to decrease with increasing T and approach 
the kinematical threshold at some critical value T = T c , 
to be identified with the deconfmement temperature. In 
this scenario one expects stable particles (poles on the 
real axis in the complex squared energy s-plane) to de- 
velop a width as a result of absorption in the thermal 
bath. At the same time, resonances (poles in the second 
Ricmann sheet in the complex s-planc) should develop 
T-dependent widths, increasing with increasing temper- 
ature. Such a resonance broadening mechanism was first 
proposed in detail in connection with dimuon production 
in heavy ion collisions [l3| . 

The natural framework to determine so has been 
that of QCD sum rules This quantum field theory 
framework is based on the operator product expansion 
(OPE) of current correlators at short distances, extended 
beyond perturbation theory, and on Cauchy's theorem 
in the complex s-planc. The latter is usually referred 
to as quark-hadron duality. Vacuum expectation 
values of quark and gluon field operators effectively 
parametrize the effects of confinement. An extension 
of this method to finite temperature was first outlined 
in [f2| . Further evidence su ppo rting the validity of 
this program was provided in llal . followed by a large 
number of applications [l6|-[l7j. Of particular interest 
to the present work are the results obtained for sq(T) 
in [TH using QCD Finite Energy Sum Rules (FESR) 
for the (light-quark) axial-vector current correlator. 
The leading dimension FESR relates sq(T) to the pion 
decay constant f v (T), and this in turn to the light- 
quark condensate (using the Gell-Mann-Oakes-Renner 
relation [H[). In the chiral limit it was found that 
s o (T)/s o (0) ~ U(T)/U(0) ~ {fMT)/(M)(0), which 
holds to a very good approximation. This relation 
hints towards the possible coincidence of the critical 
temperatures for deconfmement and chiral-symmetry 
restoration. In this paper we extend this analysis to 
finite density, thus obtaining sq{T,^,b) from FESR using 
as input the light-quark condensate at finite tempera- 
ture and density determined in the Schwinger-Dyson 
equations (SDE) framework. 

The paper is organized as follows: In Sec.|TI]we find the 
relation between the quark condensate and the PQCD 
threshold sq using FESR for the axial- vector current cor- 



relator. In Sec. IIIII we compute the quark condensate at 
finite T and fj,B from a convenient parametrization of the 
quark propagator in the SDE framework. In Sec. IIVI we 
present our analysis of the QCD phase diagram and show 
that the deconfmement and chiral symmetry restoration 
transitions take place at basically the same temperature 
to some 3 % accuracy, i.e. within the numerical precision 
of the method. We finally summarize and discuss our 
results in Sec. [Vj 



II. FINITE ENERGY QCD SUM RULES 

We begin by considering the (charged) axial-vector 
current correlator at T = 




<0|T(A M (a;)A„(0))|0), 



= -g^Il^q 2 ) + q^qMq 2 ) , (1) 

where A^(x) =: u{x)^^^,d{x) : is the axial-vector cur- 
rent, q^ = (u>,q) is the four-momentum transfer, and 
the functions 110,1 (q 2 ) are free of kinematical singulari- 
ties. Concentrating on the function TLo(q 2 ) and writing 
the OPE beyond perturbation theory in QCD [l4| . one 
of the two pillars of the sum rule method, one has 

rio(<? 2 )| Q cD = c Q i+ c 2N (q 2 , f i 2 )(d 2N ( f i 2 )) , (2) 

N=l 

where fi 2 is a renormalization scale, the Wilson coeffi- 
cients Cn depend on the Lorentz indices and quantum 
numbers of the currents and on the local gauge invari- 
ant operators On built from the quark and gluon fields 
in the QCD Lagrangian. These operators are ordered 
by increasing dimensionality and the Wilson coefficients, 
calculable in PQCD, fall off by corresponding powers 
of — q 2 . The unit operator above has dimension d = 
and CqI stands for the purely perturbative contribution. 
Hence, this OPE factorizes short distance physics, en- 
capsulated in the Wilson coefficients, and long distance 
effects parametrized by the vacuum condensates. The 
second pillar of the QCD sum rule technique is Cauchy's 
theorem in the complex squared energy s-plane 

- f ° dsf(s)ImU (s) = -^-.i dsf(s)U (s) , 
71 Jo lm Jc(\ so \) 

(3) 

where /(s) is an arbitrary analytic function, and the ra- 
dius of the circle so is large enough for QCD and the 
OPE to be used on the circle (see Fig. I). The inte- 
gral along the real s-axis involves the hadronic spec- 
tral function. This equation is the mathematical state- 
ment of what is usually referred to as quark-hadron du- 
ality. Using the OPE, Eq.@, and an integration kernel 
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FIG. 1: FESR integration contour C(|so|) in the complex 
square energy s-plane. The QCD threshold so in the FESR 
is the radius of the circle. 



f(s) = s N (N = 1, 2, ■ ■ • ) one obtains the FESR 

C 2N {6 2N )=Ak 2 f° dss N - x -lmR Q {s) 



£o_ 

N 



[l + 0{a s )] (JV = 1,2,...). (4) 



For N = 1 the dimension d = 2 term in the OPE 
does not involve any condensate, as it is not possible to 
construct a gauge invariant operator of such a dimension 
from the quark and gluon fields. Nevertheless, it is 
a-priori conceivable to generate a d = 2 term in some 
dynamical fashion, e.g. in PQCD at very high order 
(renormalons). However, there is no evidence for such 
a term (at T = 0) from FESR analyses of experimental 
data on e + e~ annihilation and r decays into hadrons 
[l9j . At very high temperatures, th oug h, there seems 
to be evidence for some d = 2 term [20|. However, the 
analysis to be reported here is performed at much lower 
values of T, so that we can safely ignore this contribution 
in the sequel. 

The extension of this pro gram to finite temperature is 
fairly straightforward [l2|, UMi with the Wilson coeffi- 
cients in the OPE, Eq.Q, remaining independent of T 
at leading order in a s , and the condensates developing a 
temperature dependence. Radiative corrections in QCD 
involve now an additional scale, i.e. the temperature, 
so that a s = a s (/i 2 , T) . This problem has not yet 
been solved successfully. Nevertheless, from the size of 
radiative corrections at T = one does not expect any 
major loss of accuracy in results from thermal FESR to 
leading order in PQCD, as long as the temperature is 
not too high, say T < 200 McV. Essentially all appli- 
cations of FESR at T ^ have been done at leading 
order in PQCD, thus implying a systematic uncertainty 
at the level of 10 %. One new feature at T ^ is 
the appearance of a new cut in the complex energy 
w-plane (l2j . and centered at the origin with extension 



— |<z| < w < \o\- This is due to a contribution to the 
current correlator in the space-like region (q 2 < 0) which 
vanishes at T = 0. Conceptually, this originates in the 
scattering of the current by either quarks (antiquarks) 
or by hadrons in the medium, in the case of QCD or the 
hadronic representation, respectively. When considering 
the rcst-framc (q — > 0) this scattering term either 
becomes a delta function of the energy or it vanishes 
identically, depending on the channel. For instance, in 
the case of the axial- vector current correlator, Eq.JT]), 
the QCD scattering term is proportional to S(lu 2 ). The 
corresponding term in the hadronic representation is 
non-zero but it is suppressed relative to the tree-level 
pion contribution, as the axial-vector current can only 
couple to an odd number of pions. Another new 
feature at finite temperature is the possible existence of 
non-scalar (Lorentz noninvariant) vacuum condensates. 
This does not affect the present analysis, as we shall 
only consider dimension d = 2 FESR. 

In the static limit (q—> 0), to leading order in PQCD, 
and for T ^ and 7^ the function U (q 2 ) in Eq.([T]) 
becomes no(u; 2 , T, /is); to simplify the notation we shall 
omit the T and /1 b dependence in the sequel. A straight- 
forward calculation of the spectral function in perturba- 
tive QCD gives 
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i Imn (s)| PQCD = —2 



1 — flA 



2 



2 



- -^T 2 S(s) [Li 2 (-e^/ T ) + Li 2 (-e-^/ T ) 

7T^ L 

where Li2(a;) is the dilogarithm function, s = ur 

1 



(5) 



and 



n±(x) 



(6) 



are the Fermi-Dirac thermal distributions for particles 
and antiparticles, respectively. We have assumed mass- 
less quarks, as quark mass corrections arc negligible. 
However, later when we determine the quark condensate 
in Section III this approximation will be relaxed. 

In the limit where T and/or hb are large, Eq. ([5]) be- 
comes 
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Imn (s)| 



PQCD 
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7T 2 T 2 



(7) 



In the hadronic sector we assume pion-pole dominance of 
the hadronic spectral function, i.e. the continuum thresh- 
old So to lie below the first radial excitation with mass 
Mjtj ~ 1300 MeV. This is a very good approximation at 
finite T, as we expect sq to be monotonically decreasing 
with increasing temperature. In this case, 
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Imn(s)| [ 



2 / 2 (7> B ) 5(.s), 



(8) 
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FIG. 2: (Color online) Lattice data and parametrization of 
the absolute value of the quark condensate as a function of T 
in the phase transition region. 



where /^(r, hb) is the pion decay constant at finite T 
and /i B , with 7^(0,0) = 92.21 ± 0.14 MeV [HJ. 

Turning to the FESR, Eq.(gJ}, with N = 1 and no 
dimension d = 2 condensate, and using Eqs.© and (JSJ 
one finds 













/ ds 


1 — n+ 


(f) 




(f)] 
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8/2(T, MB )-8T 2 [Li 2 (-e^/ T ) +Li 2 (-e-^/ T )](9) 

This is a transcendental equation determining sq(T, /j B ) 
in terms of (T, fi B ) ■ The latter is related to the light- 
quark condensate through the Gell-Mann-Oakes-Renner 
relation UM 



(10) 



/ 2 (0,0) (^)(0,0) 

where the quark and pion masses have been assumed in- 
dependent of T and fx B [22j ■ A good closed form approx- 
imation to the FESR, Eq.©, for large T and/or fi B is 

obtained using Eq.(|7|! with h + f^J ~ n_ (^-^ ~ 0, in 

which case 



a (T, ^b)^8tt 2 f n {T, MB ) - - tt 2 T 2 - 4 M | . 
Using Eq. fTU)) this can be rewritten as 

s (T,fi B ) (W>)(T,hb) (T 2 /3 - mI/tt 2 ) 



so (0,0) (#)(0,0) 



2/ 2 (0,0) 



(11) 



(12) 



The quark condensate can be computed from the in- 
medium quark propagator, whose non-perturbative prop- 
erties can be obtained e.g. from known solutions to the 
Schwinger-Dyson equations (SDE) as discussed in the 
next section. 
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mi (beV ) 


r, 


1 


-0.490 


-0.112 


2 


0.495 


0.352 


3 


-0.879 


0.259 



TABLE I: Parameters m; and ri, i = 1,2, 3 to describe the 
Lorentz covariant part of the quark propagator. 



III. QUARK PROPAGATOR AND 
CONDENSATE 



The quark condensate can be computed from the quark 
propagator S(ko,k) in Euclidean space. At finite T and 
fiB the condensate is given by 



x r d 3 k 



(13) 



(^)(T,/i B ) 
x Tr S[(2n+l)irT + i f i B ,k] 

x S[ko - (2n+ 1)ttT- ifi B ]- 
Introducing the Poisson summation formula 

5^(-l)'exp{(ifco + A*B)i/r} = 
; 

(2tt) T ^ 5[k - (2n + 1)ttT - i(x B ] 

n 

leads to 

T E j ^^S[(2n+l)rrT + lf ,B,k} = 

E(-!) ( / (0Tr[5(fc o , fc)] exp{(zfc + l* B )l/T} 

(15) 

Using this result in Eq. (fl"3]) gives 

m(T,n B ) = -N c j2(-iy J 



(14) 



x Tr[S(k , k)} exp{{iq + n B )l/T}. 



(16) 



Notice that from Eq. (jTB")) the vacuum contribution to 
the condensate comes from the term with I = 0. For this 
we use the value 



(^)| = -(0.241 GeV) 3 



(17) 
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FIG. 3: (Color online) Heat capacity for the quark condensate 
as a function of T for /is = (solid line) and fiB = 300 MeV 
(dash line). The critical temperature T c corresponds to the 
maximum of the heat capacity for a given value of \ib- 



FIG. 4: (Color online) Heat capacity for the PQCD thresh- 
old so as a function of T for /is = (solid line) and 
fiB = 300 MeV (dash line). The critical temperature T c cor- 
responds to the maximum of the heat capacity for a given 
value of fiB- 



The true matter contribution to the condensate is thus 
(#)(T,„ fl )=-iV c ^(-l) i | -0^ 



x Tr[S(k , k)] exp{(iq + n B )l/T}. 



(18) 



Due to the loss of Lorentz covariance at finite T and/or 
fiB, the general structure of the propagator is given by 



S-^ko, k) = Aj Q k + Bj-k + C, 



(19) 



where A, B and C are scalar functions of fco and k. They 
can be obtained from non-perturbative methods such as 
solutions to SDE. We adopt this procedure here. Mo- 
tivated by the success of the rainbow-ladder truncation 
of the SDE and the effective interaction of Ref. 23| in 



the description of light pseudo-scalar and vector mesons, 
and the meromorphic representation of the quark propa- 
gator p3 [. we consider the parametrization 



ijoko + ibj ■ k + 



(20) 



and choose 6, the masses mj, and the residues j-j, i = 
1 . . . 4, to be real numbers. In addition wc seek In- 
dependent values of b, and r^. The Lorentz covari- 
ant part of this parametrization is fitted by requiring 
the propagator to reproduce key features of the rainbow- 
ladder model [23| at T = 0. In particular, to match 
the ultra-violet behavior of the gap equation for massive 
u/d quarks, the value of the condensate in vacuum, and 
the constituent quark masses, as dictated by the of so- 
lutions to SDE. Table I shows the values thus obtained 



for the parameters mj and 7"j, i = 1 . . .3. The last term 
in Eq. (|20p is added to reproduce the Lorentz covari- 
ance breaking effects of the heat bath at T ^ and/or 
/ib ^ 0. The values of 6, and are adjusted to 
reproduce the light-quark condensate as a function of T 
for hb = [see Eq. (|2ip below] extracted from lattice 
QCD i by means of a point-distance minimization pro- 
cedure. Carrying out the integrations in Eq. (|18[) , and 
in terms of the parametrization of the quark propagator 
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FIG. 5: (Color online) Transition temperatures for the quark 
condensate and the PQCD threshold so as functions of the 
baryon chemical potential. 
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in Eq. (f20f , we obtain 



(V^>(7> s )| 




where K\(x) is a Bessel function, and for convenience 
we have defined hi = 1 for i = 1,2,3, and 64 = b. 
Figure [2] shows the lattice QCD data for the light quark 
condensate as a function of T Q together with the curve 
obtained from the absolute value of the sum of Eqs. (fl7|) 
and (|21[) for fis — 0. This parametrization gives a good 
description of the condensate for the range of tempera- 
tures where the phase transition occurs. 

IV. QCD PHASE DIAGRAM 

With the parametrization of lattice data at finite T 
and hb — 0, we proceed to extend the analysis to finite 
[iB- To explore the QCD phase diagram we make use of 
the expressions for the light-quark condensate and of the 
PQCD threshold sq that describe the chiral and decon- 
finement phase transitions, respectively. Next, we com- 
pute the corresponding susceptibilities which are propor- 
tional to the heat capacities, —d(tpip)/dT and —dso/dT. 
For a given fis, the transition temperature is identified 
as the value T c where the heat capacity reaches a maxi- 
mum. Figure [5] shows the transition temperatures for the 
condensate and for sq. These temperatures are basically 
identical within a small window of roughly 3 MeV around 
T = 185 MeV, for all values of /is up to the maximum 
value of hb = 300 MeV. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the QCD phase diagram 
at T 7^ and \ib 7^ based on the behavior of the light- 



quark condensate and of the PQCD threshold as probes 
of chiral symmetry restoration and deconfinement, 
respectively. We have shown that these quantities are 
related through a QCD FESR and found that they lead 
to essentially equal transition temperatures. The quark 
condensate, and thus the PQCD threshold, is computed 
using the quark propagator in the SDE framework. 
We have found it convenient to use a mcromorphic 
parametrization of this propagator in terms of real poles 
and residues. These are fixed by demanding consistency 
with the rainbow-ladder truncation of SDE at T = 0, 
and a good description of lattice QCD data for the 
quark condensate at finite T. With this simple scenario 
we have been able to extend the analysis up to baryon 
chemical potential fig ~ 300 MeV. From our results 
we can estimate the position of the critical end point 
to be hbc > 300 MeV and T c < 185 MeV, respectively. 
A more precise location of the critical end point would 
require a more refined treatment of the parametrization 
of the quark propagator. 



VI. ACKNOWLEDGMENTS 



We thank E. Cuautle and J. C. Arteaga for use- 
ful hints to carry out the parametrization procedure, 
W. Bietcnholz for pointing out useful lattice results, 
P.C. Tandy for his guidance on how to parametrize the 
quark propagator and extract quantities such as the 
condensate from this parametrization and E. Swanson 
for useful suggestions. Support for this work has been 
received in part by CONACyT (Mexico) under grant 
numbers 82230, 128534, CIC-UMSNH 4.10 and 4.22 
and PAPIIT-UNAM IN103811-3, FONDECYT 1095217 
(Chile), Proyecto Anillos ACT119 (Chile), and NRF 
(South Africa). 



[1] L. McLerran and R. D Pisarski, Nucl. Phys. A 796, 83 
(2007). 

[2] L. D. McLerran and B. Svetitsky, Phys. Lett. 98B, 195 
(1981). 

[3] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 73, 074512 
(2006); M. Cheng et al, Phys. Rev. D 74, 054507 (2006); 
Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. 
Lett. B 643, 46 (2006); Y. Aoki, S. Borsanyi, S. Durr, 
Z. Fodor, S. D. Katz, S. Krieg and K. K. Szabo, JHEP 
0906, 088 (2009); A. Bazavov et al, Phys. Rev. D 80, 
014504 (2009); S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. 
Katz, S. Krieg, C. Ratti and K. K. Szabo, JHEP 1009, 
073 (2010), S. Borsanyi et al, JHEP 1011, 077 (2010). 

[4] A. Bazavov et al, Phys. Rev. D 80, 014504 (2009); M. 



Cheng et al, Phys. Rev. D 81, 054504 (2010). 
[5] K. Kanaya, in PoS (Lattice 2010), 012 (2010). 
[6] See for example: A. Alexandru, M. Faber, I. Horvath and 

K.-F. Liu, Phys. Rev. D 72, 114513 (2005); G. Endrodi, 

Z. Fodor, S. D. Katz and K. K. Szabo, JHEP 1104, 001 

(2011). 

[7] F. Karsch et al, Nucl. Phys. Proc. Suppl. 129, 614 
(2004). 

[8] See for example G. A. Contrera, M. Orsaria, and N. 

N. Scoccola, Phys. Rev. D 82, 054026 (2010); C. A. de 

Sousa, P. Costa, M. C. Ruivo and See for example: S-H. 

Hansen, AIP Conf. Proc. 1343, 592 (2011). 
[9] See for example: S.-x. Qin, L. Chang, H. Chen, Y.-x. Liu 

and C. D. Roberts, Phys. Rev. Lett. 106, 172301 (2011); 



7 



P. K. Srivastava, S. K. Tiwari and C. P. Singh, 
larXiv:1101.lT5T] F. Xu, T. K. Mukherjee, H. Chen and 
M. Huang, EPJ Web of Conf. 13, 02004 (2011). 

[10] O. Philipsen, Nucl. Phys. A830, 713c (2009). 

[11] J. I. Kapusta and E. S. Bowman, PoS C POD2009, 018 
(2009). 

[12] A. I. Bochkarev and M. E. Shaposnikov, Nucl. Phys. B 

268, 220 (1986). 
[13] C. A. Dominguez and M. Loewe, Z. Phys. C 49, 423 

(1991). 

[14] For a review see e.g. P. Colangelo and A. Khodjamirian, 
in At the Frontier of Particle Physics/Handbook of QCD, 
edited by M. Shifman (World Scientific, Singapore 2001), 
Vol.3, p. 1495. 

[15] C.A. Dominguez and M. Loewe, Physical Review D 52, 
3143 (1995). 

[16] C. A. Dominguez and M. Loewe, Phys. Lett. B 233, 201 

(1989) . The (near) equality of the critical temperatures 
for chiral-symmetry restoration and deconfinement was 
shown analytically in A. Barducci, R. Casalbuoni, S. De 
Curtis, R. Gatto and G. Pettini, Phys. Lett. B 244, 311 

(1990) . These authors used a result for the thermal quark 
condensate valid for < T < T c , while the first reference 
only made use of the low-T expansion of chiral pertur- 
bation theory, obtaining somewhat different critical tem- 



peratures. 

[17] For recent applications see e.g. C.A. Dominguez, M. 

Loewe and J.C. Rojas, J. High Energy Phys. 0708, 040 

(2007); C. A. Dominguez, M. Loewe, J. C. Rojas, and 

Y. Zhang, Phys. Rev. D 81, 014007 (2010); ibid. D 83, 

034033 (2011), and references therein 
[18] C.A. Dominguez, M. S. Fetea, and M. Loewe, Phys. Lett. 

B 387, 151 (1996). 
[19] C. A. Dominguez and K. Schilcher, Phys. Rev. D 61, 

114020 (2000); J. High Energy Phys. 0701, 093 (2007). 
[20] E. Megias, E. R. Arriola, and L. L. Salcedo, Phys.Rev. 

D 81, 096009 (2010). 
[21] K. Nakamura et at, Particle Data Group, J. Phys. G 37, 

075021 (2010). 

[22] T. Altherr and D. Seibert, Phys. Rev. C 49, 1684 (1994). 
[23] P. Maris and CD. Roberts, Phys. Rev. C56 3369 (1997); 

P. Maris and P.C. Tandy, Phys. Rev. C60 055214 (1999); 

P. Maris and CD. Roberts Int. J. Mod. Phys. E 12 297 

(2003). 

[24] M.S. Bhagwat, M.A. Pichowsky and P.C. Tandy, Phys. 
Rev. D67 054019 (2003); R. Alkofer. W. Detmold, C. S. 
Fischer and P. Maris, Phys. Rev. D70 014014 (2004); N. 
Souchlas, J. Phys. G37 115001 (2010). 



